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To greatly enhance output of nuclear fusion produced neutrons in a laser-initiated Coulomb explo¬ 
sion of Deuterium clusters, we propose to accelerate the resulting ions by a quasi-dc electrical pulse 
to the energies where the D + + D collision cross-section is the highest. With D + ions bombarding 
then a Deuterium-rich solid-state cathode, this allows one to solve a few problems simultaneously by 
(a) completely removing electron cloud hindering the Coulomb explosion of ionic core, (b) utilizing 
up to 100% of the cluster ions to collide with the high-density packed nuclei, and (c) reaching highly 
increased cross-section of neutron production in a single D + + D collision, in particular by using a 
multi-layered target. We also consider the use of .E-pulse acceleration for diagnostic purposes. 

PACS numbers: 36.40.Gk, 52.50.Jm, 25.45.-z, 79.77.+g 


Nuclear fusion reactions in solid-deuterium laser- 
produced plasma were first observed 45 years ago [1]. A 
promising recent development in the field is generation 
of neutrons in laser-induced ionization and explosion of 
deuterium clusters [2] with a conversion efficiency up to 
~ 2 x 10 6 neutrons/J [3]. A major mechanism here is a 
Coulomb explosion, CE [2-5] of the clusters, whereby an 
irradiated and highly ionized cluster loses its free elec¬ 
trons that ideally are almost instantly swept away by the 
laser, with the ionic core torn apart by repulsive Coulomb 
forces resulting in CE. Part of the process is the forma¬ 
tion of shock-shells in expanding ionic cloud predicted 
in [6], explored in detail in [7] and most recently experi¬ 
mentally observed in [8]. In general, however, there could 
be some other mechanisms of cluster explosion, such as 
e. g. the quasi-neutral micro-plasma and hydrodynamic 
models [9], etc., more characteristic for lower laser inten¬ 
sities ( < 10 15 VE/cm 2 ). While the generation of neutrons 
indicating nuclear fusion reaction have been successfully 
observed [1-3], the results are still far from the goal of 
decades-long quest for an elusive efficient nuclear fusion 
for energy-producing application, whereby the energy of 
generated neutrons exceeds that of the input. 

Further advance in using laser-induced explosions of 
clusters is blocked by three major problems: insufficient 
kinetic energy of ions (typically a few KeV , with colli- 
sional cross-section, 07, very low), vs ~ lOOA'eE where 
Of goes up by orders of magnitude, see below; low utiliza¬ 
tion of produced ions due to low density of surrounding 
plasma (typically < 10 18 cm -3 vs ~ 10 23 cm -3 in solid 
state), hence low number of fusion collisions; and finally, 
free electrons that eventually neutralize the ion cloud and 
hamper Coulomb explosion. Shock-shells in CE [6,8] may 
increase that collision rate [6] by having ions collide in 
the higher-density area near original cluster, but this en¬ 
hancement is still insufficient to attain good efficiency. 

In this Letter we propose: (1) the way of overcom¬ 
ing those problems by using a laser-synchronized electri¬ 
cal pulse to accelerate laser-produced ions D + to ener¬ 


gies sufficiently high (up to ~ 70 — lOOA'eE) to greatly 
increase the cross-section for the neutron production, 
and smashing them against a deuterium (or tritium)-rich 
solid-state target (cathode); aside from greatly enhanc¬ 
ing neutron output, this would provide an ultimate test 
of laser-l-cluster nuclear fusion for energy production ap¬ 
plications, (2) the enhancement of the solid-state target 
performance by making it as a stack of thin layers, and 
(3) a diagnostic of the cluster and explosion structure for 
research purposes and application as a neutron source. 

This approach may be viewed as a cross between laser- 
induced CE and basic electrostatic generation of neutrons 
[12] providing substantial cross benefits. First, there is no 
need anymore to strive for powerful laser irradiation, in 
particular to remove the electron cloud; the laser energy 
has to be just enough to attain a reasonable ionization, 
not to produce high ion energies or blow electron cloud 
away. Ionized electrons here are quickly removed from 
the expanding cloud, thus barring them from neutraliz¬ 
ing the ions in the cloud and hindering the useful effect 
of CE. Similarly, laser-induced ionization may allow to 
use lower E-field than in electrostatic generators. The 
system also simplifies the analysis and diagnostics of the 
process: once electrons are removed, the dynamics of re¬ 
maining cloud of positive ions is strictly due to a repulsive 
Coulomb explosion, which now may differ from an ideal 
CE only in that in each case the radial density distribu¬ 
tion of the ions may be non-uniform depending on their 
initial distribution. The latter one can thus be elicited by 
using segmentation of the cathode into a few electrically 
isolated sections/rings and recording the time-dependant 
current from each one of them, see below. 

In a common arrangement, a cluster or a jet of clusters 
is injected between an anode collecting ionization pro¬ 
duced electrons and a cathode covered by deuterium or 
tritium-rich material, as a target for D + ions accelerated 
by a strong A-pulse applied to the electrodes. The ion en¬ 
ergy then can get much higher than that produced by CE 
and reach the optimum domain of up to 100AeE; it is di- 
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rectly controlled by the impulse amplitude. A solid-state 
target insures then a high probability of fusion collisions. 
An .E-pulse has to be sufficiently long to be maintained 
till all the ions reach the cathode. Its duration for e. 
g. 50KeV voltage and 2 cm electrode spacing, has to be 
> 10ns, while the laser pulse is typically sub -ps long; the 
formation time of ion cloud is even much shorter. For 
a laser intensity 10 15 bF/c?n 2 , free electrons are pulled 
from the core faster than a laser cycle [6], and then they 
are swept away by the E-pulse and brought to the anode. 

Electrode geometry may vary from spheres to cylin¬ 
ders, and to cones, while a parallel configuration provides 
for the simplest arrangement, Fig.l. Consider a cluster 



FIG. 1: Full ion flow rate, ( dN/dr ) / TVs vs dimensionless 
time, r = t/to, for various potential/kinetic energy ratios, 
X = Uo/Tomax ■ Inset: the schematics of the proposed experi¬ 
ment; dotted circles outline the edge of ion cloud after cluster 
explosion at two different moments. The expanding positively 
charged cloud moves down toward the deuterium-rich cathode 
and finally hit it under the action of high-voltage potential. 


where Ny is a total number of ions in a cluster, rije is the 
ion charge (for hydrogen atoms or isotopes, rij = 1), and 
R 0 is an original radius of cluster. When R ci » R 0 , 
where R c i(t) is a cloud radius, the ion motion is unaf¬ 
fected by ion collisions, and the ion movement is iner¬ 
tial in the frame of the center of mass (COM), while 
COM accelerates toward a cathode with its 2 -speed be¬ 
ing v COM (t) = —Uot/zoM, with t = 0 at the moment of 
explosion. Similarly to the Hubble expansion, the ion ra¬ 
dial distance p from COM is proportional to their original 
velocity, p = v 0 t, v 0 < v 0ma x = \/2T 0 max /M. 

By far, the most crucial “boost-factor” in the proposed 
scheme is the ion acceleration due to quasi-dc poten¬ 
tial. We consider the simplest neutron-producing reac¬ 
tion due to collision of two deuterium nuclei: d + d —> 
3 He(0.82MeV) + n(2A5MeV). The importance of suf¬ 
ficiently high nuclei energy transpires from the fact that 
the fusion cross-section cr/(ej n ), increases with ei n by 
many orders of magnitude. Here ej n = Co + Tq is a 
full energy of the collision of fast ion D + with D target 
molecules (such as e. g. E 2 , CD 4 , or D 2 0 [10,11]; D 2 
is the most efficient one [14]). Following [11] and using 
fusion probability defined as 


w f (e in ) 



NtarVCTf {e)de/ ( de/dt ), 


( 2 ) 


where de/dt is the energy loss rate largely due to nucleus- 
nucleus scattering (the contribution of neutron genera¬ 
tion to this rate could be ignored), N tar is the number 
density of target nuclei, and v is the velocity of deuterons. 
In the range of energies of interest and regular solid-state 
density, Wf is approximated by a simple formula [11] 

~ Cnexp f* \Je sc / (3) 


initially at the distance zq from a cathode and a potential 
between them is Co = gVo, where Vo is a respective volt¬ 
age and q - an ionic charge (for E + , q = e). Once bound 
electrons are freed by laser pulse, the E-pulse “vacuum¬ 
cleaning” pulls them to the anode. The ensuing ionic CE 
would then typically produce a shock at the edges of a 
CE cloud [6-8], yet soon the density of over-run ions de¬ 
creases. Then, for an ideal CE, the expansion proceeds as 
if a sphere with uniform density distribution in the mo¬ 
mentum space. In general, however, this could be more 
complicated, due to expansion originated by thermal ex¬ 
plosion, secondary ionization, or compositions, such as 
e. g. non-uniform or heterogeneous clusters [3,13] of 
different ionic species, or mixed clusters formed by de¬ 
positing layers of different atoms, etc. This may result 
in distinctly different initial non-uniform ion kinetic en¬ 
ergy radial distribution, with the same maximum energy 
Tomax , see below. In the case of an ideal CE, whereby 
the initial density is almost uniform, we have [6] 

Tomax = {eni ) 2 Ny/Rq, (1) 


where Cd = 0.18 and e sc = IMeV. Eq. (3) is best fitted 
for D 2 molecules; considering an example of cluster with 
Ny ~ 3,250, and Rq = a(3IVi;/47r) 1 / 3 with a ~ 10 _8 cto 
being an initial averaged spacing between ions (assuming 
an ideal case whereby all the atoms ionized by laser), we 
have T 0max ~ 5 KeV, Eq. (1). In such a case, if one 
chooses Co ~ 100 KeV, Wf is boosted by ~ 0.5 x 10 13 . 
This makes it clear how far a Coulomb explosion is from 
producing a substantial neutron output. However, as im¬ 
pressive as this number may look, the probability w / per 
se is not the bottom-line to strive for, as far as energy nu¬ 
clear fusion is concerned. For that, the minimum goal is 
to have an averaged output energy of neutrons, eoWf, to 
exceed that used to ionize deuterium and accelerate D + 
ions, 6 i n , so that the “in/out” efficiency 77 exceeds unity, 
77 = eoWf/ein > 1, where en = 2.45MeV is the energy of 
a neutron exiting fusion reaction. For real applications 
the input power “at wall-plug” must be higher, but 77 > 1 
provides a minimal requirement. The efficiency 77 vs in¬ 
put energy ei n is depicted in Fig. 2; one can see that 77 
never reaches unity; it peaks at e; ra = e sc /4 = 1.75MeH 









3 


[i. e. actually outside the domain of validity of Eq. (3)], 
and even then rjmax = ^(eoCo /e sc )e~ 2 ~ 0.033 <C 1. For 
e sc < 200-ftfeE, we have 7 ] < 0.006, and for e sc < 100 KeV 
< 10 -3 . Thus, the efficient fusion is unreachable here; 
essentially, this is true for any beam-target fusion, in¬ 
cluding laser-accelerated ion beams [15]. 



FIG. 2: The “in/out” efficiency of neutron production for 
bulk, r), and stack (multi-layered), r/st, targets vs incident 
energy of deuteron nucleus. Inset: a schematic configuration 
of the stack target; e st is the total potential across the stack. 


Barring the target compression (as in inertial con¬ 
finement fusion, with a required number density up to 
N ~ 10 25 cm -3 [10]), one way to enhance the target per¬ 
formance is to have it as a stack of thin multi-layers to 
be crossed by ion beam, and apply a moderate dc voltage 
between them to replenish the energy lost by ions. This 
proposed “cascading” or “recycling” is to mitigate the 
rapid ion energy loss in a bulk target. This loss results 
in dramatic reduction of fusion cross-section cr/(e) in Eq. 
(3), so that a good part of the bulk’s volume is lost for 
the neutron production. Making the layers thin enough 
to have a low energy loss Ae <C ej„ per layer (to be re¬ 
plenished by an equal inter-layer potential), and their 
number respectively large, N st = e st /Ae 1, where e et 
is the total potential across the stack, we have the prob¬ 
ability of fusion per layer as (dwf /dei n ) Ae, and the total 
cascade probability w s t ~ e s t(dw // dei n ), such that 


r) st = 


tpWst 

£st T Gn 


rj{ein) 



tst/2 


tin f-st + G: 


( 4 ) 


If we choose e. g. e st = Cj n , the cascade target in the 
realistic setting offers the enhancement of \Je sc /ei n / 4. 
over a bulk target. In that case, r] s t(ein) peaks at €i n = 
e sc /9 « 778 KeV, with {j] st )max = 0.023; r] st = 0.0021 at 
£in = lOOiveE, and 0.0086 at €i n = 200 KeV. In spite 
of enhancement, we still have r) st < 1, i. e. the efficient 
fusion is still unreachable. Similar calculations for a D + 
ions and Tritium target, (with even higher output neu¬ 
tron energy, ~ 14.1M eV, instead of 2A5Mev for D + +D) 
leave this conclusion stand. As was mentioned, the same 
is true for any beam+target fusion mode, as opposed to 


the inertial confinement fusion, whereby the efficiency 
can in principle exceed unity. The major reason for this is 
that the product of density and confinement time in the 
former case is far below the one required by the so called 
Lawson criterion [10,11]; there is simply not enough time 
for an D + ion to enter into a ?r-producing reaction in a 
bulk target [16], while in a cascade target it needs energy 
replenishment; the situation here is remotely reminiscent 
of a critical mass phenomenon in fission reactions. A po¬ 
tentially viable path for energy solution could be a hybrid 
approach combining ion-beam and inertial confinement 
fusion, but it is outside the scope of this Letter. A pos¬ 
sible way to go might be bombarding the solid-state tar¬ 
get by small ionized clusters, instead of separate ions, or 
chunks of partly-decomposed clusters sufficiently accel¬ 
erated by a E-pulse, which may create longer-sustained 
high-temperature at the impact area (see also below), 
which we are planning to address elsewhere. 

Regardless of energy-producing prospects, the pro¬ 
posed system may be used both as a portable source of 
neutrons, with the advantage of having orders of mag¬ 
nitude better time-clocking of the output neutron pulse 
over electrostatic sources [12] due to well synchronized 
laser pulse trigger, and a greatly useful research tool 
for diagnostics of cluster structure and explosion cloud. 
While the diagnostics of ultra-short pulse laser-produced 
“macro”-plasma have been extensively studied (see e.g. 
[17]), the nanoscale mapping of fs pulse absorption was 
studied [18] recently using a ’’plasma explosion imag¬ 
ing” . Yet combined with E-pulse, the system may def¬ 
initely offer new possibilities. Its diagnostics capability 
comes naturally from the fact that it has already the ma¬ 
jor components (CE + E-field) employed by a so called 
Coulomb Explosion Imaging [19-22] (CEI), used first in 
1989 [19] to yield images of small individual molecules. 
The CEI technique has been developed into a fine and 
sophisticated tool [20], including most recent use of pixel¬ 
imaging mass-spectrometer camera [21], and attaining 
the first image of the Efimov trirner in helium [22]. How¬ 
ever, there is an important difference with cluster CE: a 
cluster is comprised of too many atoms/ions, and aside 
perhaps from studying the details of its surface, which 
might be of interest to the physics of cluster formation, 
one is dealing with much more “macro-effects” than those 
in the small individual molecules imaging, and might be 
less concerned about imaging. As an example, an impor¬ 
tant subject to explore is the cluster composition, to be 
elicited from analyzing the time dynamics of total ionic 
current at the cathode, Fig. 1, as well as ’’differential” 
currents flowing through sub-cathodes, [23]. The advan¬ 
tage also is that usually a single laser shot and single 
cluster as e. g. in [8] and [18] can be used instead of 
averaging over many shots as in CEL 

To have an idea of expected characteristics for those 
applications, we consider some details of ion flow dy¬ 
namics. We introduce a potential/kinetic energy ratio, 
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X = Uo/To max , and a dimensionless time, r = t/to, 
where to = zq \/2M/Uq is a time for COM to reach a 
cathode due to potential Uq alone. The first ions hit the 
cathode in tim e T m j n a nd the last ones - in T max , where 
Tmin,max = a/ 1 + X -1 T At/2; At = 2/y/x and At is 
total duration of ion flow. Ions with the same starting 
energy To make an expanding sphere (for its outer edge 
see inset, Fig. 1); which falls down to a cathode with 
acceleration —Uq/zoM. The rate number of ions hitting 
the cathode at T m in < t < Tmax is as: 


1 dN 
Ns dr 


^ (f + 72 ) F W; F = /JVm, (5) 


where £ = To/To max is a relative initial kinetic energies of 
ions, and function /(£) describes a radial distribution of 
these energies that satisfies a condition /q/(£)^(£ 3 ^ 2 ) = 
1. In the case of ideal CE, we have fcE = const = 1. 
To illustrate density/energy-profile sensitivity of this sys¬ 
tem, we will also consider two other distinctly differ¬ 
ent models of that distribution, in particular “hot ball”, 
/hb = 5£/3, which has a hollow core, while its outer shell 
is populated by hottest ions, thus making it a sustained 
shock, and a “cool ball”, fcB = 5(1 — £)/2, with a dense 
cold core and a vanishing density of “hot” outer ions. 

In Eq. (5), the integration limits £i(r) and £ 2 ( 1 ") are 
determined by the area of the cathode engaged. If all the 
ions hitting cathode are included, we have £2 = 1 , and 
Ci = Cmin(i") = (x/4) (t -1 - t) <1, which is a minimal 
initial energy of ions reaching a cathode at the moment r, 
so that for CE, Fce = 1 — £mn, for hot ball, Fhb = 5(1 — 
Cmin)/6, and for a cool ball, Fcb = 5(1 - £mn) 2 /4. For 
the CE case, the rate dN/dr in units Ns for x from 2.5 
to 25 is depicted in Fig. 1 (if T 0max = AKeV, that would 
correspond to Uq ranging from 10 to 100 KeV). While 
initial kinetic energy of an ion is To, which increases to 
T = U 0 + T 0 when it hits a cathode, the total energy of 
the cloud delivered to the cathode during entire process 
in an ideal CE case, is Ts = Ns(Uo + 3T 0max /5). 

To prevent ions from hitting an anode, one needs a suf¬ 
ficient overhead spacing, z up , between clusters and anode, 
(with total cathode-anode spacing zca = zo + z up > z cr , 
and thus a sufficient voltage between the plates, Vca > 
V cr ), so that Zcr/zo = V cr /Vo = 1 + X • The max¬ 
imum “hot spot” radius p sp at the cathode is reached 
at T sp , which are respectively as: p sp /zo = 2^1 + x/x> 
T sp = 2y -1 + 1, i. e. the spot gets tighter as x i n_ 

creases, as expected. The effect is sensitive to the ion 
energies distribution /(£) in the explosion, and thus may 
offer well-resolved time-of-flight diagnostics. It can be 
implemented by segmenting the cathode into isolated 
concentric rings or sub-cathodes, and recording ion flow 
in each of them, as well as a total count of the ions, for 
various distributions /(C); for the example of 3 rings, see 
[23]. Ultimately, the cathode can be made as a large set 
of tiny pixels (with resolution determined then by pixel 


numbers), as well as to serve as a mass-spectrometer tool 
in search and production of multi-atomic ” nano-chunks” 
to be explored in the attempt to bridge ion-beam and 
inertial confinement fusion. 

In conclusion, we proposed to reach ultimately high 
yield of neutrons in laser-driven explosion of deuterium 
clusters, by using a synchronized T-pulse with up to 
100 KeV peak field, and making D + ions bombard a neg¬ 
atively charged deuteriunr-rich cathode. The maximum 
yield is reached as free electrons are removed from the 
ion cloud, D + ions are fully utilized and made to collide 
with solid-state deuterium-rich target instead of plasma. 
We also proposed to further enhance the output by using 
a multi-layer target. While the energy production goal 
appears to be still unreachable, the major application 
of the system could be an efficient neutron source with 
laser-controlled timing of neutron pulses. The system can 
also be made into a sensitive diagnostic tool to resolve an 
intrinsic structure of ion cloud and cluster itself. 


* Electronic address: alexander.kaplan@jhu.edu 

[1] F. Floux, D. Cognard, L.-G. Denoeud, G. Piar, D. 
Parisot, J. L. Bobin, F. Delobeau, and C. Fauquignon, 
Phys. Rev. A 1 , 821 (1970). 

[2] T. Ditmire, J. Zweiback, V. P. Yanovsky, T. E. Cowan, G. 
Hays, and K. B. Wharton, Nature (London) 398, 489492 
(1999); J. Zweiback, R. A. Smith, T. E. Cowan, G. Hays, 
K. B. Wharton, V. P. Yanovsky, and T. Ditmire, Phys. 
Rev. Lett. 84, 2634 (2000); G. Grillon, Pli. Balcou, J.- 
P. Chambaret, D. Hulin, J. Martino, S. Moustaizis, L. 
Notebaert, M. Pittman, Th. Pussieux, A. Rousse, J.-Ph. 
Rousseau, S. Sebban, O. Sublemontier, and M. Schmidt, 
Phys. Rev. Lett. 89, 065005 (2002). 

[3] H. Y. Lu, J. S. Liu, C. Wang, W. T. Wang, Z. L. Zhou, 
A. H. Deng, C. Q. Xia, Y. Xu, X. M. Lu, Y. H. Jiang, Y. 
X. Leng, X. Y. Liang, G. Q. Ni, R. X. Li, and Z. Z. Xu, 
Phys. Rev. A80, 051201 (2009). 

[4] Gotts, N. G., Lethbridge, P. G. and Stace, A. J. J. Chem. 
Phys. 96, 408-421 (1992); K. W. Madison, P. K. Patel, 
M. Allen, D. Price, T. Ditmire, J. Opt. Soc. Am. B 20, 
113 (2003). 

[5] M. Lezius, S. Dobosz, D. Normand, and M. Schmidt, 
Phys. Rev. Lett. 80, 261 (1998); G. A. Mourou, T. 
Tajima, and S. V. Bulanov, Rev. Mod. Phys., 78, 309 
(2006); Th. Fennel, K.-H. Meiwes-Broer, J. Tigges- 
bumker P.-G. Reinhard, P. M. Dinh and E. Suraud, ibid,, 
82, 1795 (2010); K. I. Popov, V. Yu. Bychenkov, W. Roz- 
mus, and L. Ramunno, Phys. Plasmas 17, 083110 (2010). 

[ 6 ] A. E. Kaplan, B. Y. Dubetsky, P. L. Shkolnikov, Phys. 
Rev. Letts, 91, 143401 (2003). 

[7] F. Peano, R. A. Fonseca, and L. O. Silva, Phys. Rev. 
Lett. 94, 033401 (2005); V. Yu. Bychenkov and V. F. 
Kovalev, Plasma Phys. Reports, 31, 178 (2005); 

[ 8 ] D. D. Hickstein, F. Dollar, J. A. Gaffney, M. E. Foord, 
G. M. Petrov, B. B. Palm, K. E. Keister, J. L. Ellis, C. Y. 
Ding, S. B. Libby, J. L. Jimenez, H. C. Kapteyn, M. M. 
Murnane, and W. Xiong, Phys. Rev. Letts, 112, 115004 
(2014). 

[9] A. V. Gurevich, L. V. Pariskaya, and L. P. Pitaeveskii, 








5 


Sov. Phys. JETP, 22, 449 (1966), also 36, 274 (1973); 
J. E. Allen and J. G. Andrews, J. Plasma Phys. 4, 187 
(1970); L. M. Wickens, J. E. Allen, and P. T. Rumsby, 
Phys. Rev. Lett. 41, 243 (1978); M. A. True, J. R. Albrit¬ 
ton, and E. A. Williams, Phys. Fluids 24, 1885 (1981); 
H. M. Milchberg, S. J. McNaught, and E. Parra, Phys. 
Rev. E64, 056402 (2001). 

[10] S. Atzeni and J. Meyer-ter-Vehn, “The Physics of Inertial 
Fusion”, Oxford University Press, 2004. 

[11] V. P. Krainov and B. M. Smirnov, Sov. Phys. JETP, 105, 
559 (2007). 

[12] J. Reijonen, K.-N. Leung, R. B. Firestone, J. A. En¬ 
glish, D. L. Perry, A. Smith, F. Gicquel, M. Sun, H. 
Koivunoro, T.-P. Lou, B. Bandong, G. Garabedian, Zs. 
Revay, L. Szentmiklosi, G. Molnar, Nucl. Instrum, and 
Meth., A522, pp.598-602, (2004); J. Csikai, Ed. “Hand¬ 
book of Fast Neutron Generators”; CRC Press: Boca 
Raton, FL, 1987; Vols. I and II, 242 and 224 pp. 

[13] H. Pauly, “Atom, Molecule and Cluster Beams” 
(Springer, Berlin, 2000), Vol. II, Sec. 2.6. 

[14] H. S. Bosch and G. M. Hale, Nucl. Fusion, 32, 611 (1992). 

[15] For a review of ion acceleration by ultra-intense lasers, see 
e. g. M. Roth, E. Brambrink, P. Audebert, M. Basko, A. 
Blazevic, R. Clarke, J. Cobble, T. E. Cowan, J. Fernan¬ 
dez, J. Fuchs, M. Hegelich, K. Ledingham, B. G. Logan, 
D. Neely, H. Ruh, and M. Schollmeier, Plasma Physics 
and Controlled Fusion 47, B841 (2005) 

[16] Be reminded that a neutron generation reaction requires 
a collision by two deuterium nuclei, which due to huge 


repulsing Coulomb potential is much less probable then 
atom-ion collisions between D and D + . 

[17] M. Roth, J. Instrumentation, 6, R09001 (2011). ICAL 
REVIEW A 91, 053424 (2015) 

[18] D. D. Hickstein, F. Dollar, J. L. Ellis, K. J. Schnitzen- 
baumer, K. E. Keister, G. M. Petrov, C. Ding, B. B. 
Palm, J. A. Gaffney, M. E. Foord, S. B. Libby, G. 
Dukovic, J. L. Jimenez, H. C. Kapteyn, M. M. Murnane, 
and W. Xiong, ACS Nano, 8, 8810 (2014). 

[19] Z. Vager, R. Naaman, E. P. Kanter, Science, 244, 424 
(1989). 

[20] H. Hasegawa, Chem.Phys.Lett. 349, 57 (2001); C. Cor- 
naggia, Laser Physics, 19, 1660 (2009); H. Xu, T. Okino, 
K. Nakai, and K. Yamanouchi in ’’Progress in Ultrafast 
Intense Laser Science VII”, K. Yamanouchi, D. Char- 
alambidis, D. Normand, Eds, Springer Series in Chemical 
Physics 100, 35 (Heidelberg, 2011); 

[21] C.S. Slater, S. Blake, M. Brouard, A. Lauer, C. Vallance, 
C. S. Bohum, L. Cristensen, J. Nielsen, M. P. Johansson, 
H. Stapelfeldt, Phys.Rev. A 91, 053424 (2015); 

[22] M. Kunitski, S. Zeller, J. Voigtsberger, A. Kalinin, L. Ph. 
H. Schmidt, M. Schoffler, A. Czasch, W. Schollkopf, R. E. 
Grisenti, T. Jahnke, D. Blume, and R. Dorner, Science, 
348, 551 (2015). 

[23] See Supplemental Material at [URL will be inserted by 
publisher] for the demonstration example of the diagnos¬ 
tics of the cluster density profile for different models of 
those profiles using the set of three sub-cathodes. 



6 


SUPPLEMENTAL MATERIAL for 

Coulomb cluster explosion boosted by a quasi-cfc pulse diagnostic tool and ultimate 

test of laser fusion efficiency in clusters 

A. E. Kaplan 

Electr. and Comp. Eng. Dept, The Johns Hopkins University, Baltimore, MD 21218 

(Dated: November 27, 2015) 


This Section is to illustrate how the diagnostic of clus¬ 
ter density profile can be attained by segmenting the 
cathode into isolated concentric rings or sub-cathodes, 
and recording ion flow in each of them, as well as a total 
count of the ions, for various distributions /(£), see the 
text in two paragraphs following Eq. (5). The calcula¬ 
tions here are based on Eq. (5) in the MAIN SECTION, 
where we consider three model radial density profiles: 

(a) an ideal Coulomb explosion, fcE = const = 1, 

(b) “hot ball” profile, Jhb = 5£/3, with a hollow core, 
while its outer shell is populated by hottest ions, thus 
making it a sustained shock, and 

(c) a “cool ball” profile, fcB = 5(1 — £)/2, with a dense 
cold core and a vanishing density of “hot” outer ions. 

For the demonstration purposes we consider here the set 
of three rings, the central one being a disc with a radius 
p < p\, a middle ring pi < p c < P2, and external ring 
P 2 < Pc < Pmax > where we set the sizes in such a way 
that in the end of the process, the total ion flow in each 
of them be the same for all three sub-cathodes for the 
ideal CE case, whereby they have to be pi/p sp ~ 0.475 
and P 2 Ip S p ~ 0.709. 

Fig. 1 depicts the time dynamics of the current/flow 
of ions through each ring and total ion count (see in¬ 
serts) in each of those rings for each chosen model. It 
illustrates a great potential of such a system for spec¬ 
troscopy and diagnostic purposes. The ion flow dynamics 
for various models in Fig. 1 clearly indicates that the ef¬ 
fect is sensitive to the distribution of ion energies /(£) in 
the original explosion, and thus may offer a well-resolved 
time-of-flight diagnostics of that distribution. Depend¬ 
ing on application, the electrodes geometry can be made 
as spherical or cylindrical surfaces. In general case, in¬ 
stead of rings, the cathod can be made of multiple pixels 
that would provide for much greater temporal and spatial 
resolution of explosion ion cloud and its dynamics. 





FIG. 1: (a) Rates of the current/flow of ions through sub¬ 
cathode rings (1 - inner, 2 - medium, 3 - outer rings), 
and time-integrated ion count (inserts) in each of those 
rings for various models: (a) an ideal Coulomb explosion 
(CE), with even distribution of ions inside the cloud; (b) 
a “hot ball”, with a hollow core and outer shell of hot 
ions, and (c) a “cool ball”, with a dense colder core and 
vanishing density of “hot” outer ions. 






